Sub-Planck scale structures in a vibrating molecule in the presence of decoherence 
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We study the effect of decoherence on the sub-Planck scale structures of the vibrational wave 
packet of a molecule. The time evolution of these wave packets is investigated under the influence 
of a photonic or phononic environment. We determine the master equation describing the reduced 
dynamics of the wave-packet and analyze the sensitivity of the sub-Planck structures against deco- 
herence in the case of a hydrogen iodide (HI) molecule. 
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I. INTRODUCTION 

Recent progress of controlled femtosecond pulses has 
advanced greatly the technology during the last few years 
[l]. A new field of molecular optics has emerged where 
lasers are used to manipulate the internal and exter- 
nal degrees of freedom of molecules, to deflect beams of 
molecules, to control molecular dynamics, and to align 
molecules 0, Q. Many investigations have focused on 
the vibrational motion of diatomic molecules. The single 
bond between the atoms acts as a spring and supports 
harmonic oscillations for small amplitudes, but the bond 
can break (dissociate) when stretched too much. These 
phenomena usually occur at time scales between few pi- 
coseconds and few hundred femtoseconds. With ultra- 
short pulses one can now prepare a molecular wave packet 
and probe its evolution and observe molecular reactions 
in this time domain. Successful experiments have been 
performed on several molecules [4] . The most convenient 
model for studying the vibrational motion of diatomic 
molecule is the Morse potential, which is an exactly solv- 
able system (Hj. Coherent superposition of several vi- 
brational levels of the molecule creates the wave packet 
which, due to quantum interference, shows revival and 
fractional revivals 0, @] in their time evolution. 

Fractional revivals are associated with superpositions 
of separated wave packets (for example, the so-called 
Schrodinger cat states), which manifest clear quantum 
interference effects and nonclassical features, which can 
be well visualized in the phase space of the vibrational 
motion. Although the experimental observation of small 
quantum interference structures is very challenging, it 
has already been visualized in Pico- meter scale [9]. A 
number of different phase space distribution functions 
have been introduced [10] and investigated over the years, 
and among these the Wigner distribution [ll|] is particu- 
larly useful, because its negativity yields an indication of 
nonclassical behavior [l2|, [ij, [l^] . Zurek [l^ first showed 
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that this negativity reveals the existence of the smallest 
structures in phase space i.e., the sub-Planck scale struc- 
tures. One may expect that Heisenberg's uncertainty 
principle implies that structures on scales smaller than 
the Planck constant have no observable consequence, 
while instead Zurek [is] showed that these highly non- 
classical structures are expected to be particularly sen- 
sitive to decoherence. Through a short walk in contro- 
versy, recently sub-Planck scale structures draw consider- 
able attention and have been found by others in different 

situations [il, [13, E H S Ell, S H, m, m, [13 . 



Decoherence due to the coupling to an external envi- 
ronment is the main responsible for the disappearance 
of nonclassical manifestations of quantum states and it 
is considered one of the mechanisms through which the 
classical world at the macroscopic level emerges from 
the quantum substrate 0, [2^, [sO] • Decoherence on the 
molecular vibrational degree of freedom is due to the cou- 
pling between vibrational and rotational modes [31] and 
also to the coupling with the photonic and phononic de- 
grees of freedom The latter are associated with a 
super- Ohmic environment describable in terms of a con- 
tinuous set of bosonic modes and in this paper we shall 
focus on the effect of this source of decoherence. To be 
more specific, we will study the sub-Planck scale struc- 
tures in the Morse system and the effect of decoherence 
on these structures in molecular wave packets. We shall 
determine the master equation describing the reduced 
dynamics of the wave packet and analyze the robustness 
of the sub-Planck structures against decoherence. 



The paper is organized as follows. In Sec. II we give 
a brief overview of the Morse potential and its coherent 
states, while in Sec. Ill we derive the master equation in 
the case of the coupling with a bosonic environment at 
thermal equilibrium. In Sec. IV we study the effect of de- 
coherence on the Wigner function at the sub-Planck level 
and the sensitivity to decoherence of these structures is 
analyzed. Finally, we conclude in Sec. V. 
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II. REVIEW OF THE MORSE MODEL OF A 
VIBRATING MOLECULE 



III. MASTER EQUATION FOR THE MORSE 
OSCILLATOR 



Vibrational dynamics of diatomic molecules are well 
described by Morse potential [HI, [s^, IH, [s^ . It can 

be described as 



V{x) = D{e 



-2(3x 



(1) 



where x = r/ro — 1, ro is the equilibrium value of the 
inter- nuclear distance r and /3 is a range parameter. D is 
the dissociation energy, which has been extensively stud- 
ied in a wider context of this model [s^, • Defining 
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and 



8^ 



(2) 



where /i is the reduced mass of the vibrational motion, 
the eigenfunctions of the Morse potential can be written 



(3) 



where ^ = 2Ae-^^, < ^ < oo, and n = 0, 1, [A - 1/2], 
with [p] denoting the integer part of p, so that the total 
number of bound states is [A — 1/2]. The parameters A 
and s satisfy the constraint condition s + 2n = 2A — 1. 

Note that A is potential dependent and s is related to 
energy E and, by definition, A > 0, 5 > 0. In Eq. (|3|), 
Lf^{y) is the associated Laguerre polynomial and N is the 
normalization constant: 



TV 
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Revival and fractional revivals appear during the time 
evolution of a suitably prepared wave packet and are well 
studied in the literature p, 0, Q • Here we study the ef- 
fect of decoherence on the motion of a molecular wave 
packet through its sub-Planck scale structures. There 
structures are found at one eighth of the fractional re- 
vival time in the Wigner phase space distribution [2^. 
The initial wave packet is taken here as SU(1, 1) coher- 
ent state (CS) of this potential [so'], which is obtained 
upon applying the displacement operator on the ground 
state. The CS is given by 



oo 
k=0 



1) 



fc!r(i 
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where < k < [X — 1/2] correspond to the bound states 
of the Morse potential and k > [A — 1/2] are the appro- 
priate scattering states. The parameter 77 is associated 
with the "amplitude" of the CS and possesses the same 
phase of the displacement amplitude a, while its mod- 
ulus is given by the relation I77I = tanh|a| ^]. In our 
numerical analysis, we will always consider low energy 
coherent states well below the dissociation limit so that 
only the bound states of the Morse potential can be used 
as basis set. 



As described in Sec. I, we now investigate the effect of 
the decoherence of an external phononic or photonic en- 
vironment on the sub-Planck scale structure. Therefore, 
the total model Hamiltonian is [2^1 



H — Hs 



Hp 



Hi, 



(6) 



where Hgys is the Morse Hamiltonian of the vibrational 
mode. He is the environment Hamiltonian described by 
a set of independent bosonic modes 



He = ^hujk{alak 



1/2), 



(7) 



and Hi is the interaction between the Morse particle and 
the environment, which we choose of the following form 
(see also 



Hi 



H.c, 



(8) 



where ak are coupling constants. This choice corresponds 
to assume the rotating wave approximation (RWA) in 
the interaction with the environment so that we neglect 
counter-rotating terms, while the operator O is a generic 
operator of the vibrational mode, whose specific form 
depends on the considered environment. Using standard 
techniques [29], one gets in the usual Born-Markovian 
approximation, the following master equation for the re- 
duced density operator of the Morse oscillator p. 



d_ 



-jXHsys.p] (9) 

OspO^ + dpd\ - d^d2P - p6\6 
6\p6 + o^pOx - 66\p - p6y6^\ , 



where the operators Oj {j = 1,2) are new operators of 
the vibrational mode corresponding to "modifications" 
of the operator associated to the absorption from the 
environment (Oi) or emission into the environment (O2) 
of vibrational quanta. This fact is easily understood if 
we look at their expression in the energy eigenbasis |n, s) 
used in Eq. ([5]). In fact, one has 

Oi=E^r"Ks)(^>5i. (10) 



O"^'"Trg{L0mn)\(r{uJmn)\^n{L0mn), (H) 



where 

= 0™■"7^^;(a;™„)|a(u;™„)^n(u;„^„) + l].(12) 

The quantities O^'^ are the matrix elements of O, 
gi^mn) is the density of states at the transition frequency 
between two energy levels, ujrnn = {Em — En)/h, and 



3 



(a) 



30 



P 



30 




(b) 



Max 



30 



P 



30 




■0.3 0.3 0.6 

X 



■0.3 0.3 0.6 

X 



Min 



FIG. 1: (Color online) Time evolution of Morse wave packet in phase space Wigner distribution: (a) cat state at time t = jTrev 
and (b) sub-Planck scale structures, appeared at the middle at time t = |Trev, where a = 0.3, f3 = 2.07932, ro = 3.04159 
a.u., and n = 4. Here, x and p are the dimensionless position and momentum variables, where x = r/ro — 1 and p is the 
corresponding scaled variable. 



^i^mn) = [expjficjmn/^B^} ~ IS the mean thermal 
number of environmental excitations, being the latter at 
equilibrium at temperature T. The appearance of these 
two operators is a direct consequence of the nonlinearity 
of the molecular vibrational motion. In fact, in the linear 
case the transition frequencies ujrnn do not depend on n 
and m and therefore Oi and O2 becomes proportional to 
O. As a consequence, master equation (j9j) become iden- 
tical to the master equation of a harmonic oscillator in a 
thermal environment in the RWA [29]. 



IV. SUB-PLANCK SCALE STRUCTURE AND 
ITS SENSITIVITY THROUGH DECOHERENCE 

We now solve the master equation (|9]) for the specific 
case of the HI molecule and we adopt the Wigner function 
picture in order to look at the effects of decoherence on 
sub-Planck scale structures in phase space. The Wigner 
distribution is defined {h = 1) by 

W{x,p,t) = ^l^ {x-^\p{t)\x+^)e"'Pdx', (13) 

and well describes the nonclassical interference effects as- 
sociated with the time evolution of a wave packet in the 
nonlinear potential of the Morse oscillator. 

This fact is visible, for example, in Fig.[Tl which shows 
the time evolution of an initial CS wave packet in phase 
space at two different fractional revival times in the ab- 
sence of decoherence. We have considered a HI molecule, 
which has 30 bound states, with (3 = 2.0793, reduced 
mass^= 1819.99 a.u., ro 3.0416 a.u., and D 0.1125 
a.u. [23]. We have assumed here (and also in the follow- 
ing) that the initial wave packet is well below the dissoci- 
ation limit, so that it involves only the lower levels of HI 
molecule (the energy distribution is peaked around the 



n = 4 vibrational level). Figure [Tf a) shows the vibra- 
tional cat state after one fourth of the fractional revival 
time. Here, the revival time is Trev = 4.89 x 10^ a.u. 
Due to the anharmonicity of the system, one can notice 
the different squeezing effects in the two separated CSs 
forming the cat state. The number of ripples in the inter- 
ference region increases for increasing mean energy of the 
initial CS. The sub-Planck scale structures appear in the 
interference region at one eighth of the fractional revival 
time [Fig. [TJb)], where one has a coherent superposition 
of four well distinct states, forming a so-called compass 
state [15]. For this reason we shall focus our attention on 
the effect of decoherence at this fractional revival time. 

In the case of a molecular vibration, a bosonic envi- 
ronment well describes either the coupling via the dipole 
interaction with the outside electromagnetic field or, in 
the case of a molecule immersed in a liquid or gas, the 
coupling with the acoustic modes of the solvent. In both 
cases the operator O is connected with the position op- 
erator of our Morse oscillator. In fact, O describes the 
upper triangular part (in the energy basis representa- 
tion) of the dipole moment operator of the molecule in 
the electromagnetic case and of the vibrational coordi- 
nate X in the acoustic phonon bath case. However the 
two situations are analogous because the dipole moment 
operator is proportional to x. Both environments are 
super- Ohmic, that is, we have 

27ra'^{u)g{u) = Sco^, (14) 

with S characterizing the strength of the system- 
environment coupling. The physical meaning of the pa- 
rameter 5 can be seen from the fact that the master equa- 
tion ([9]) implies that the relaxation rate from level i to 
level j, Tij at zero temperature is given by 

ri,=544, (15) 
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FIG. 2: (Color online) Wigner function of the coherent state at one eighth of revival time, approximately equal to a compass 
state, for three different values of the decoherence parameter S (a) 5 = 0; (b) 6 = 0.54 x 10^ a.u.; and (c) (5 = 2.2 x 10^ a.u. 
Here, x and p are the dimensionless position and momentum variables, where x = r/ro — 1 and p is the corresponding scaled 
variable. The environment temperature is fixed at T = lOhuJoi/kB- 



where r^j is the corresponding matrix element of the po- 
sition operator between the two vibrational levels. Here 
we have chosen the coupling constant 6 such that the ra- 
tio Foi/cJoi ranges from 0.5 x 10"^ to 12.5 x 10"^. These 
values correspond to reasonable values of the coupling 
constant S; in fact, considering a typical electric dipole 
moment of a diatomic molecule one gets S :^ 10~^^ m~^ 
= 4.78 a.u. for the case of an electromagnetic envi- 
ronment. Instead, considering a dilute solvent one gets 
^ 2 X 10~^^ m~^ s^ = 95.66 a.u. for the case of a 
phononic environment. One should note that position 
variable in our study is a dimensionless quantity (scaled 
by ro). Hence, in the case of an electromagnetic environ- 
ment, the order of 6 in our case would be ^ = 4.37 x 10^ 
a.u., which is consistent with our study (see Fig. [3]). 
Moreover, the temperature of the environment is kept 
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FIG. 3: Comparative variations between the left (dotted line) 
and right (solid line) peaks and the central negative sub- 
Planck region (dashed line) at ^Trev with the coupling pa- 
rameter 6 (in unit of 10^ a.u.). The environment temperature 
is T = lOhuJoi/kB. 

fixed at T = lOhcooi/kB- Figure [2] shows the Wigner dis- 
tribution at one eighth of the fractional revival time for 
different values of the coupling with the bosonic environ- 
ment. Figure [2l^a) refers to no decoherence {5 = 0) and 
therefore corresponds to Fig. [T](b). Figure [2l^b) instead 



corresponds to S = 0.54 x 10^ a.u. and Fig. [2](c) corre- 
sponds to a stronger decoherence, S = 2.2 x 10^ a.u. One 
can clearly see that by increasing the coupling with the 
bosonic environment, the interference region is more and 
more affected. 
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FIG. 4: Variation in the central negative sub-Planck region 
at ^Trev with the coupling parameter S (in unit of 10^ a.u.). 
Dots are the numerical data from our analysis. It satisfies 
an exponential law (solid line) Ae~^^ , with A = 0.5847 and 
c = 0.3585. The environment temperature is T = lOhoooi/kB- 

As for the harmonic oscillator case [l5|, , decoher- 
ence affects the structure as a whole; also here for the 
Morse oscillator, the sub-Planck structures due to quan- 
tum interference are more affected than the individual 
isolated coherent state components. In fact, a distinct 
difference can be observed in the decay rate of the am- 
plitude of the sub-Planck scale structures and of the in- 
dividual CSs. This is quantitatively shown in Fig. [3l 
where these decay rates are plotted versus the decoher- 
ence strength S. We consider the left and right peaks 
of the CSs at p = and a negative peak at x = 0.077 
andp = —6.064, appearing in the sub-Planck interference 
region. The plot shows that the sub-Planck scale struc- 
ture, i.e., the central interference patterns (dashed line 
in Fig. [3|), disappears earlier compared to the individual 
CSs, as it happens in the harmonic case. It is possible 
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FIG. 5: (Color online) Variation in the central negative sub-Planck region at time ^Trev with the environment temperature T 
in the case S = 0.54 x 10^ a.u. It follows a Bose-distribution law, aexp{-6/[e^"/^ - 1]}, for a = 0.5799 and b = 0.0127. Inset 
of the figure shows the variation near the critical temperature (Tc = 0.6688). 
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FIG. 6: (Color online) Wigner distribution at times (a) ^Trev and (b) ^Trev, with coupling parameter 6 — 0.72 x 10^ a.u. and 
temperature T = IO/icjoi/^b- Central interference region in (a) (sub-Planck scale structure) disappears at a later time in (b). 
Here, x and p are the dimensionless position and momentum variables, where x = r/ro — 1 and p is the corresponding scaled 
variable. 



to see that the decay of the ampHtude of the sub-Planck 
structure follows very well an exponential law as a func- 
tion of the decoherence strength (5, as expected in usual 
bosonic environments [28]. We find that a linear expo- 
nential function Ae~'^^ well fits with our results, with 
A = 0.5847 and c = 0.3585. Figure H shows how the 
rate of amplitude damping of the chosen negative inter- 
ference region (i.e., sub-Planck region) matches with the 
exponential form. 

It is now worth seeing the effect of environment tem- 
perature on decoherence for a fixed value of the coupling 
constant 5. Owing to Eqs. ([TT]) and ([12]), one expects a 
Bose-Einstein dependence on temperature of the decay 
of the interference structures associated with sub-Planck 
structures, aexp{— 6/[e^^/^ — 1]}? where Tc corresponds 
to an effective transition temperature below which the 
discrete structure of the energy levels of the Morse os- 
cillator starts to manifests itself. This is confirmed by 
Fig. [5l where the numerical results for the value of the 
negative peak are plotted versus temperature. The data, 
corresponding to ^ = 0.54 x 10^ a.u., are well fitted by 
the above curve and the optimal fitting parameters are 
a = 0.5799, b = 0.0127 and = 0.6688. The data follow 



an exponential decay for T/Tc ^ 1, while the devia- 
tion from the exponential law (associated with the Bose- 
Einstein distribution dependence) is clearly visible only 
at very low temperatures, T < Tc, in the magnified view 
in the inset of Fig. [5l 

So far we have been studying the decoherence effect 
on the sub-Planck scale structures at 1/8 fractional re- 
vival time. Hence, it is a natural question to ask what 
happens at larger times when one can also obtain sub- 
Planck scale structures in the interference region of four- 
way break up of a coherent state. Thus, we extend 
our study to the four way break-up or the decoherence 
through sub-Planck scale regions at 3/8 and 5/8 frac- 
tional revival times. One expects a larger infiuence of 
decoherence on the sub-Planck structures for increasing 
times and this is confirmed by Fig.[6l Interference fringes 
in phase space are still visible at 3/8 fractional revival 
time, while in Fig. [BJb), corresponding to 5/8 fractional 
revival time, one can see that the sub-Planck structures 
completely disappear due to the larger decoherence ef- 
fect, whereas the individual coherent states remain al- 
most intact. The environment temperature is kept con- 
stant at T = lOhuJoi/kB, the coupling constant being 
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5 = 0.72 X 10^ a.u. 



V. CONCLUSIONS 

We have investigated the time evolution of a coherent 
state wave packet in the Morse potential under the influ- 
ence of a bosonic environment describing either photonic 
or phononic excitations. We have studied the effect of 
decoherence on the sub-Planck structures in phase space 
by looking at the evolution of the Wigner distribution. 
As it happens for the harmonic case, sub-Planck scale 



structures come out as the most sensitive to decoherence. 
A quantitative analysis provides an exponential decay of 
the amplitude of the quantum interference structures as a 
function of the coupling with the environment, in agree- 
ment with usual predictions [28]. Influence of the envi- 
ronment temperature on the decoherence is also shown 
quantitatively. This is according to the Bose-distribution 
law. Longer time effect on the decoherence is shown for 
providing another way to see the sensitiveness of sub- 
Planck scale structures compare to their original coun- 
terparts. 
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